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Abstract: Topological techniques are used to study the motions of systems of 
point vortices in the infinite plane, in singly-periodic arrays, and in doubly-periodic 
0^ ■ lattices. The reduction of each system using its symmetries is described in detail. Re- 

stricting to three vortices with zero net circulation, each reduced system is described 
by a one degree of freedom Hamiltonian. The phase portrait of this reduced system 
I is subdivided into regimes using the separatrix motions, and a braid representing the 

1—5 ' topology of all vortex motions in each regime is computed. This braid also describes 

Q\ ■ the isotopy class of the advection homeomorphism induced by the vortex motion. 

^ ! The Thurston-Nielsen theory is then used to analyse these isotopy classes, and in 

^ I certain cases strong conclusions about the dynamics of the advection can be made. 

> '■ 

^ . ^1 Introduction. 

Q , The modelling of incompressible flow at high Reynolds number as a potential flow with 

I embedded vortices has repeatedly proven useful both for analytical and numerical purposes. 
^ ■ The subject has inspired numerous reviews, each stressing different aspects of the field. 
0\ '. The articles by Aref [Afl], Chorin [Cl-3], Leonard [L], Majda [Mj], Moffatt & Tsinober 
i [MfT], PuUin & Saffman [PS], Saffman & Baker [SB], Saffman [S], Sarpkaya [Srp], Shariff & 

^1 Leonard [SL] , and Zabusky [Z] provide a representative sample. We are concerned here with 
the further simplification of modelling a two-dimensional flow by a finite collection of point 
^ . vortices. While this system is admittedly highly idealised, it has found application and, to 
r-| I some extent, experimental verification starting with von Karman's analysis in 1912 of the 
J-^ I instability of the vortex street wake behind a cylinder and Onsager's 1949 explanation of the 
^ ■ emergence of large coherent vortices in two-dimensional flow through an 'inverse cascade' of 

^ , energy. (See the literature cited for further details.) 

^ i In this paper we study the evolution of three kinds of point vortex systems, in the infinite 

■ ■ ' plane, in singly-periodic arrays, and in doubly-periodic lattices. An array of point vortices is 
useful in modelling flows in shear layers, wakes and jets. The motivation for studying vortex 
lattices comes from the problem of two-dimensional turbulence, a paradigm of atmospheric 
and oceanographic flows, and to a lesser extent from the study of vortex patterns formed in 
superfluid Helium. 

Kirchhoff recognised that the evolution of N point vortices could be formulated within the 
framework of classical mechanics as an A-degree of freedom Hamiltonian system. In many 
ways point vortices are the fluid mechanical analog of point masses evolving under the mutual 
interaction of Newtonian gravity. As with the A-body problem, point vortex systems have 
continuous symmetries that give rise to integrals of motion. In the case of zero net circulation 
(a case of considerable physical interest), the two integrals arising from translation invariance 
are in involution. Thus for three vortices, the Hamiltonian system is completely integrable. 
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Fixing a value of these integrals and factoring out by the corresponding group action gives rise 
to a new Hamiltonian system with two less degrees of freedom. For three vortex systems, this 
process of Jacobi-Poincare reduction yields a single degree of freedom system. Vortex arrays 
and lattices also have a discrete symmetry coming from the array or lattice structure. For 
rational circulation ratios these result in a discrete symmetry of the reduced system whose 
phase space is consequently a cylinder or torus. Symmetries and reduction are explained in 
detail in §3. 

Since the reduction is accomplished by factoring out the translational invariance, tra- 
jectories in the reduced system describe the evolution of the shape and orientation of the 
triangle spanned by the three vortices. Because the reduced system has one degree of free- 
dom the generic bounded orbit is periodic. This implies that the corresponding triangle 
of vortices resumes its shape and orientation after one period. However, the triangle may 
be translated. This translation vector is a dynamic phase of the type that arises when the 
evolution of a full system is reconstructed from a periodic trajectory of a reduced system. 
(For a general introduction to reduction and dynamic phases see [MZl] and [MZ2]. Further 
references are given in §3.) 

Having performed the reduction on a three- vortex system, the dynamics of the vortices 
can be understood by analysing a one degree of freedom Hamiltonian system (see [Af2], 
[ASl] and [SA]). In such systems the phase portrait is organised by the saddle points and 
the connections between them, known as separatriccs (sec Figures 3.1 and 3.2). For vortices 
the reduced Hamiltonian will also have singularities or poles arising from collisions, or more 
properly, from the superposition of pairs of vortices as collisions do not occur when there 
is zero total circulation. The separatriccs naturally divide the phase space into regions 
(called regimes in this paper) in which one would expect the corresponding fluid motions 
to share certain dynamical characteristics. These characteristics are topological and are 
precisely described using Artin's braid group. The assignment of a braid to a vortex motion 
is accomplished by visualising the vortex motion in a three dimensional space-time with 
the two space directions in the usual xy-planc and time going upwards. The motion in the 
vortex frame is periodic, and so the strands representing the motions in this frame have the 
same initial and final position. The translation into a mathematical braid in Artin's group 
is accomplished by projecting onto a fixed plane which is orthogonal to the xy -plane. This 
construction is described in §4 where it is shown that all the vortex motions arising from 
the same regime are described by the same braid. Examples in the infinite plane and in an 
array are worked out in detail. 

As a Hamiltonian system the vortex evolution takes place in a 2n-dimensional phase 
space. However, from a fluid mechanical perspective it is much more natural to view the 
motion as N trajectories in the plane (or the cylinder or torus in the array and lattice 
cases). Since the vortices carry concentrated vorticity, they generate a velocity field in the 
surrounding fluid. The evolution of passive tracer particles in this velocity field is called 
advection. In mathematical terms, advection is the evolution according to solutions of the 
differential equation corresponding to the vector field of velocities. Returning to the case of 
three vortices with zero net circulation, we pass to the vortex frame to eliminate the dynamic 
phase and so obtain a periodic motion of the vortices. In this frame the vortices generate 
a periodic velocity field, and so we define the advection homeomorphism as the Poincare 
map obtained by advecting for one period in this frame. Thus iterates of the advection 
homeomorphism describe the dynamics of advection in the vortex frame (see §6.1). 
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In §5.1 the braid of a regime is connected to the corresponding advection homeomor- 
phisms using the notion of isotopy. Two homeomorphisms are isotopic if one can be obtained 
from the other by a continuous deformation. In addition to its role as a description of the 
topology of a periodic motion of points in the plane, a braid on N strands also naturally 
describes an isotopy class of homeomorphisms on the A'^-punctured plane. Since all the mo- 
tions in a regime have the same braid, they all generate isotopic advection homeomorphisms. 
This knowledge is then applied using Thurston-Nielsen theory of surface homeomorphisms. 

The Thurston-Nielsen theory (described in §5.2 and §5.3) contains a classification the- 
orem for isotopy classes of surface homeomorphisms. In the case of a pseudoAnosov class 
the theory gives dynamical behaviour that must be present in every homeomorphism in the 
isotopy class. This information includes a lower bound for such quantities and structures 
as the growth rate of the number of periodic orbit points as the period grows, the topo- 
logical entropy, the growth rate of the length of topologically nontrivial material lines, and 
the topology of invariant manifold templates. In particular, any homeomorphism in a pA 
class will be chaotic under any of the usual definitions of the word. For three vortices in the 
plane with zero net circulation, the advection homeomorphisms are never in pseudoAnosov 
classes, but for vortex arrays and lattices, pseudoAnosov classes are common. Examples are 
discussed in detail in §6.2 and §6.3. 

Because we believe the methods described in this paper have general applicability, an 
effort has been made to make the paper accessible to a general scientific reader, perhaps 
lacking expertise in topology and/or fiuid mechanics. A fair amount of expository material 
has been included. In certain cases (eg. §4.1 and §5.1) basic material needed in a form not 
available in the literature is described in some mathematical detail. The reader may prefer 
to skip these sections in a first reading. 

§2 Equations of motion and Hamiltonians. 

We shall be considering the evolution of three types of systems of vortices: finite collec- 
tions in the infinite plane, singly-periodic arrays and doubly-periodic lattices. The names 
"array" and "lattice" will be reserved for the singly- and doubly-periodic systems, respec- 
tively. Both real and complex notation will be used for points in the plane with z — x + iy 
in all cases. A singly-periodic array of N vortices of period (or width) L is a collection of 
vortices in the plane with exactly vortices in each vertical strip of width L and if there is a 
vortex at position given by a complex number Za there is also one at Za + nL for all integers 
n. The definition of a vortex lattice requires a pair of linearly independent complex num- 
bers LVi and u!2 (the half-periods) that determine the double-periodicity. A doubly-periodic 
lattice of vortices with lattice structure generated by ui and U2 is a collection of point 
vortices that has exactly A'^ vortices in each fundamental parallelogram determined by 2uJi 
and 2u!2, and if there is a vortex at position Za there is also one at Za + ni2uji + n22uj2 for 
all integers ni and n2- Eventually we shall describe vortex arrays or lattices by singly- or 
doubly-periodic coordinates, or equivalently, as N vortices on a cylinder or torus. However, 
initially the arrays or lattices are specified by a collection of A^ distinguished points. Thus, 
for example, {zi + L,Z2, zs) represents a different array than {zi, Z2, zs). Note also that our 
definition of lattice requires that it maintain the same periodicity for all time, so a uniformly 
rotating system is not a lattice in our sense. 

§2.1 Equations of motion for the vortex systems. The a*^ point vortex has a constant 
circulation given by the real number Fq, and a position at time t given by the complex number 
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Za(t). Using an asterisk denotes the complex conjugate, the N differential equations that 
describe the evolution of the vortex systems are 

f -^Er.^.fe-.,) (2.1) 

for a = 1, ... ,N. Using A* = for the infinite plane, A" = ^ for arrays, and A" = £ for 
lattices, the complex- valued function (px in the various cases is 

<Pp{z) = - 
z 

Mz)-l^ot{^z) (2.2) 

where in the lattice case, ({z) = ({z; 001,002) is the Weierstrass zeta function, A is the area of 
a fundamental parallegram, and S = ^ — ^7^- Note that in each case the function is odd 
and has the required periodicity; (f)£{z + riL) = (j)c{z) and 4>a{z + 2niuJi + 2n2U}2) = 4'a{z) 
(the latter identity uses the Legendre relation for the Weierstrass zeta function) . 

In the presence of the evolving vortex system, a passive particle with position given by 
z{t) advects according to 

with the function the same as above. One may think of the passive particle as a vortex 
with zero circulation, and so the advection problem is the analog of the Newtonian (A^ + 1)- 
body problem. 

The equations of motion for the plane and array case are well known. For a derivation 
of the equations for vortex lattices see [O] and [SA]. The function (j) is usually thought of as 
representing the contribution of a single vortex. It has a simple pole at the position of the 
vortex and thus contributes a delta function to the curl of the velocity field. Note, however, 
that in the lattice case the function (f)j^ is not meromorphic. This is a refiection of a property 
of doubly-periodic meromorphic (i.e. elliptic) functions, namely, the residues of the poles in 
a fundamental parallelogram must sum to zero (equivalently, by Greens' Theorem, a doubly- 
periodic vector field must have zero net curl in each fundamental parallelogram). Thus the 
concept of a single doubly-periodic vortex or of a lattice (as we have defined it) with non-zero 
net circulation requires additional consideration. On the other hand, if X] Tq = 0, as will be 
assumed below, the right hand side of (2.3) is an elliptic function of z and one may speak 
without contradiction of the cfi'cct of the entire lattice on a passive particle. 

§2.2 The Hamiltonian framework. It is well known that the equations of motion for 
systems of point vortices can be put into the framework of classical Hamiltonian mechanics. 
For an A?"- vortex system the phase space is — T;^ where T;t' is the collision set defined in 
the various cases as 

T-p(2;) = {(-^1, ■ ■ ■ , zjsi) e C"^ : Zi — Zj for some i ^ j} 

'^Aiz) — {(-^1) ■ ■ ■ ) ^n) £ C'^ : Zi — Zj + nL for some i ^ j,n & Z} 

^jc{z) — {(-^1) ■ ■ ■ ) ^n) G C"^ ■ Zi — Zj + 2nui + 2rmjj2 for some i j,n,m & Z}. 
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The X and y positions of each vortex, and t/a, are conjugate variables, so the dynamical 
system is 2A^-dimensional. The Hamiltonian system has degrees of freedom, but there is 
no configuration space in the usual sense. Put in geometric language, the Hamiltonian is 
defined on a 2A'"-dimensional symplectic manifold that is not a cotangent bundle. 
In all cases the real- valued Hamiltonian takes the form 

Hx{zi, ... ,ZN)^-^Y.'^aTp^x{Za-Zp) (2.4) 

where the sum is over all a and and the primed summation symbol indicates that the 
case a = 13 is excluded. The function (essentially the real part of the antiderivative of 
(px) in the various cases is 



^v{z) =log(|^| 



,7r 



^A{z)=\og{\siYi{-z)\) ^2.5) 

<^c{z)^\og{\a{z)\) + Re{^)-^zz^ 

with a{z) the Weierstrass sigma function with half periods uJi and UJ2- The equations of 
motion (2.1) then have the form 

at dya OXa 

The presence of the factor Fq, on the left hand side indicates that in order to put the 
equations in the proper Hamiltonian form we need to change coordinates or else use a 
slightly nonstandard symplectic form. We choose the latter course and adopt the form 

^TadxaAdya, (2.7) 
and so the skew-symmetric matrix representing the form is 



-r 
r 



with r = diag(rQ,). Thus if W — {xi, . . . , xn, yi, ... , yjv, ), the equations of motion (2.2) 
become 

dt 

with V the real gradient. 

It is clear from equations (2.1) that the complex quantity J = J2^aZa is a constant 
of the motion. This integral is associated with the invariance of the Hamiltonian under 
simultaneous planar translation of all the vortices. Writing J — Q + iP with P and Q 
real- valued, and using the symplectic form (2.7), the Poisson bracket of P and Q is {P, Q} = 
(VP)"^ (fi"^)-^ VQ = J2^a- Throughout most of this paper we will be assuming that 
J2^a = 0, and so P and Q are in involution. The gradients of P and Q are obviously 
independent. Thus in the case N — 3 the system (2.6) is completely integrable. The 
structure of the level sets is discussed in §3.5 below. 



5 



The advection equation (2.3) may also be put in Hamiltonian form. Assuming that a 
motion {zp{t)} of the vortices is given, the time-dependent real- valued Hamiltonian is 

1 ^ 

Gx{z) = E r/3$;^(^ - Zfs{t)). (2.8) 

/3=1 

In this case we use the standard symplectic form on the plane. The function Gx is more 
commonly called the stream function of the advection problem. The right hand side of (2.8) 
is the velocity field of an unsteady two-dimensional fluid motion generated by the evolving 
collection of vortices. 

One can also combine the systems (2.4) and (2.8) into a single (A^-|- 1) -degree of freedom 
system, but we do not pursue this point of view here. 

§3 Symmetry and reduction. 

The evolution of the vortex systems is certainly independent of the choice of the origin. 
This is reflected in the invariance of the Hamiltonian (2.4) under simultaneous translation of 
all the vortices. There is a well-known process going back to Jacobi and Poincare which uses 
such symmetries to reduce the effective dimension of a Hamiltonian system. This process 
of reduction proceeds by first fixing a value of the integrals arising from the symmetries and 
then identifying elements in the level set that correspond under the restricted symmetry. 
The resulting reduced system is Hamiltonian and has one less degree of freedom for every 
independent continuous symmetry. The mathematical theory of reduction has been much 
developed in recent years (see, for example, [M], Appendix 5 in [A2], or Chapter V.D in 
[MH]). We shall need only a small piece of the theory here. Reduction for vortices in the 
infinite plane is described in greater generality in [AR]. 

In addition to the continuous symmetries, vortex arrays and lattices have discrete sym- 
metries that describe their periodicities. The utilisation of the discrete symmetries is rather 
different than continuous ones, resulting in a change of the topology of the phase space rather 
than a reduction of dimension. 

§3.1 Group actions. Symmetries are usually described by invariance under a group action. 
If G is a group (which we write additively), an action of G on a space X is a (at least 
continuous) map G x X ^ X, where g ■ x is written for the image of {g,x). The action is 
required to satisfy ■ x = x and {gi + g2) ■ x = gi ■ {g2 ■ x). The orbit or trajectory through x 
is the image of x under all the group elements, sometimes written G ■ x = {g ■ x : g E G}. 

Given a Hamiltonian H : X ^ R and an action of the group G, the Hamiltonian is 
said to be invariant under the group action if H{g ■ x) = H{x) for all x & X and g & G. 
Because of this invariance the Hamiltonian will be constant on the orbits of the action. Thus 
if we create a new space by identifying all the points on the same orbit, we will get a new 
Hamiltonian function on the quotient space. This process of collapsing orbits to points is 
called factoring or mod-ing out by the group. 

A group action is also generated by the collection of solutions to an autonomous ordinary 
differential equation, x = F{x) with a; G X (if both F and X are sufficiently well behaved). 
This action of M on X is usually written ip : X x R X and it satisfies dilj{x,t)/dt = 
F{ip{x, t)). In Dynamical Systems this action is called a flow, and here we use fiow exclusively 
in this sense. The evolution of a fluid is called a fluid motion. Note that if a fluid motion 
is steady in a region X, then the evolution of passive scalars (advection) is described by a 
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mathematical flow. However, if the velocity field F depends on t, the motion of the fiuid is 
not a fiow in the sense used here. One can obtain a flow on a new space X x Rhy defining 
a new variable r G M with f = 1 and letting x = F{x, r). 

For Hamiltonian systems, invariance under a group action is closely related to integrals of 
motion: if H is invariant under an action of M and this action is the same as the flow induced 
by another Hamiltonian if, then H is an integral invariant of the original Hamiltonian 
system. If there are k independent symmetries that commute (usually expressed as an 
action of r'^), then there will be k independent integrals of motion. 

§3.2 Continuous and discrete symmetries of the vortex systems. The Hamiltonian 
(2.4) is invariant under simultaneous translation of all the vortices, 

Hx{zi+T, ... ,ZN + r) = H{zi, ... ,zn) (3.1) 

for T e C. This action of C (or K^) is called the continuous symmetry. It gives rise to the 
integral of motion J := Q + iP deflned in §2.2. 

In addition to the continuous symmetries, vortex arrays and lattices have a symmetry 
associated with their periodicity. For arrays, 

Hj^{zi+niL, ... ,ZN + nNL) ^ H^{zi, ... ,zn) (3.2) 

for all collections of integers ni,n2, ... ,nN. Thus the Hamiltonian is invariant under an 
action of Z^. In the lattice case the Hamiltonian is invariant under a (Z^)^ action, 

Hjr{zi + ni2uji + mi2uj2, ... , zn + nN^tui + mN'2uj2) = Hc{zi, ... ,zn)- (3.3) 

We shall call these the discrete symmetries. 

§3.3 Reduction and reconstruction using formulae. From this point onward we restrict 
to the case of three vortices, = 3. In addition, the sum of the circulations is assumed to 
be zero, I] T/J = 0. 

li Z — zi — Z2, then using the definition of J and the fact that 5] Fq = 0, 

-J + FiZ 

Z2- - 



Z\ - Z2. 



-./ %,Z P.4) 

Substituting this into the result of subtracting the first two equations of (2.1) yields 

f .Z|(,,(,),,,(Z™),,,(^)). (3.) 

This is a one degree of freedom system with Hamiltonian 

K,(Z) = Il(*^(Z) + £^t,(Z£^) + ^it^(^i^^^)) (3.6) 

using the standard symplectic form (phase portraits in a planar and array example are 
shown in Figures 3.1 and 3.2, respectively). Thus the two independent integrals Q and P 
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Figure 3.1: Phase portrait of the reduced Hamiltonian system for 
three vortices in the infinite plane with circulations 1, 1/2 and —3/2. 
The upper case letters label regimes discussed in §4.3. 




Figure 3.2: Phase portrait of the reduced Hamiltonian system for 
three vortices in a singly-periodic array with circulations 1, 1/2 and 
—3/2. The Roman numerals label regimes discussed in §4.5. 



have reduced the system by two degrees of freedom, from three to one. However, the way 
in which the continuous symmetries have brought about the reduction is not clear. One 
would expect that Kx could be obtained by substituting (3.4) into the Hamiltonian (2.4), 
but there is an unexplained factor of — ^^jf^. 

When Ps is an integer the discrete symmetries in the array case are reflected in the 
invariance of Kj[ under translations by P3L. Thus we may think of the reduced system as 
living on a cylinder of width P3L. But although zi and Z2 have a width L symmetry, why 
should their difference have a width F^L symmetry? These questions are answered using the 
geometric framework of the next subsection. 
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Reconstruction refers to the process of obtaining solutions to the full equations (2.1) 
from the reduced one (3.5). After fixing a value of J, this may be done by "quadrature". 
Explicitly, using (3.4) in the differential equation (2.1) for zi, 



dzi 



2m 



1 



(12 cPx{Z{t))+Vz <Px{ 



-J - TaZ 



))■ 



(3.7) 



Thus once an initial position 2;i(0) is chosen, zi{t) may be obtained by integrating (3.7), and 
then Z2{t) = zi{t) - Z(t) and (3.4) yields z^it). 

Note that if Z{t) is given and the initial position of the first vortex -2i(0) is changed by a 
translation to 2:1 (0) = z\{Qi) + T for r e C, then the entire motion translates by r, i.e. the new 
solution are Za{t) = Za{t) + r, for a = 1,2, 3. Also note that if Z{t) is periodic with period 
P, (3.7) shows that the vortex positions Za{t) will, in general, not be periodic. Rather there 
will be a constant fe e C with Zct{t + P) = Za{t) + b, for all t and a. The number b is an 
example of a dynamic phase which is also best understood in the geometric framework of 
the next subsection. 

§3.4 Reduction in the geometric framework. For most of the remainder of this section 
we restrict to the case of singly-periodic arrays and drop the subscripts in various notations, 
so (f) — (f)j[, etc. For expositional simplicity, let L — 1. 

A natural starting point is the observation that the system does not change if we pick 
an integer n and change Za to Za + n. This means that we may treat each Za as a singly- 
periodic variable, i.e. an element of a cylinder x R. This is equivalent to factoring out 
by the (z)"^ symmetry and obtaining a Hamiltonian system with phase space {S^ x m)"^ — T. 
Unfortunately the integral Q — ^a^a is not well defined as a function of the singly-periodic 
x-coordinates. Put another way, the continuous symmetry of the ]R^-action is still present in 
the factored system, but it is no longer Hamiltonian in character; the one-parameter subgroup 
that translates in the |/-direction is not a Hamiltonian fiow on the product of cylinders. Thus 
we proceed by first reducing the system by the continuous symmetries and then using the 
discrete symmetries. 

Fix a value of the integrals, say Jo = Qq + IPq, and let Aq be the integral manifold (level 

set) 



It is convenient to use coordinates {Z, Y) for Aq with Z = ai — a2 and Y = ai for a = 
(01,02,03) e Aq. Since J2^a{(ia + t) = J2^a da, the set Aq is invariant under the action 
of 'M? as well as under the Hamiltonian fiow. In coordinates, r e acts as {Z, Y + t). We 
now mod out by the M? action restricted to Aq, i.e. we consider a and a -|- t(1, 1, 1) to be 
the same point; in coordinates this just amounts to identifying Y with Y + r for all r G C. 
The reduced space, denoted R, is thus a copy of the complex plane coordinatized by Z. 
The Meyer-Marsden-Weinstein Theorem (see [M]) assures us that R has a unique symplectic 
form that pulls back under the projection — > it! to the restriction of of the form Q, on 
Aq. To compute this form on R, note that restricted to Aq, J2^adxa = = ^I'adya- Then 
using E Tq = 0, 



Aq = {{zi,Z2,Z3,) -.^TaZ, 



Jo}. 



^adxa A dy, 



-rir2 



{dxi — dx2) A {dyi — dy2) 



-riFa 



dZ^ A dZ'^ 
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where Z = + iZ'^ . Thus the proper symplectic form on R is times the standard 

one. The flow restricted to projects to a Hamihonian flow on /?, and the Hamihonian is 
obtained by substituting (3.4) into the Hamihonian (2.4) of the fuU system. As noted above, 
this substitution yields (3.6) without a factor — ^-p^, a factor now explained by the proper 
symplectic structure on R. 

The next step is to connect the Hamiltonian flow on R with the full flow on (this is 
the reconstruction step). The flow on may be expressed as 

(Z(i;Zo),y(i;yo,^o)) (3.8) 

with the symbols after the semicolon representing the initial conditions. Note that the 
evolution of the Z coordinate just depends on the Z initial condition and not on Y . The 
invariance of the Hamiltonian under the M^-action is expressed by 

{Z{t- Zo), Y{t; Yo + r, Zq)) = {Z{t; Zq), Y{t; Yq, Zq)) + r. (3.9) 

In more concrete terms, the variable Z = zi — Z2 in conjunction with the flxed value 
of J determines the shape and orientation of the triangle determined by the vortices. The 
variable Y = zi gives its position. Translating this position by r has only the effect of 
translating the entire evolution by r. 

Since the Hamiltonian system on R has just one degree of freedom, the generic bounded 
orbit is periodic. Thus after some period P, the shape and orientation of the triangle of 
vortices will be reestablished, but in general, the entire conflguration could have translated 
by some amount. This amount is the dynamic phase and is equal to 

b{Zo) := Y{P; Yo, Zq) - F(0; Fq, Zq) (3.10) 

with Zq an element of the chosen periodic orbit. By virtue of (3.9) this expression is in- 
dependent of Yq. In geometric language this phase represents the return map (holonomy) 
in the flber of the bundle — > -R as one goes around the periodic orbit loop in R. The 
numerical value of the phase may be obtained from Z{t) by integrating (3.7). 

Thus far we have focused on the role of the continuous symmetries in reducing the 
number of degrees of freedom of the vortex systems. The use of the discrete symmetries in 
the array and lattice systems will be called factoring to distinguish it from reduction. For 
expositional simplicity we now assume that Fi = 1 and will freely use this numerical value as 
well as the symbol in various formulas. This results in no loss of generality because equations 
(2.1) show that rescaling all circulations by the same amount only changes the speed (and 
perhaps the direction) of the vortices, but does not change their trajectories. 

The main observation required in using the discrete symmetries in the reduction process is 
that additional symmetries can only be used if they preserve the level sets Aq, or equivalently, 
the value of the integral J . If (ai, ^s) is in Aq there is no guarantee that (ai + ni, 02 + 
^2,03 + n^) is in Aq for arbitrary (ni, 712,^3). However, it is easy to flnd conditions that 
insure this, namely, T^^ani = 0. The set of all such (ni, 712,^3) is a subgroup of 1? which 
we denote K. This subgroup always contains all multiples of (1, 1, 1) and is larger exactly 
when Fs is rational. Since translation by integer multiples of (1,1,1) is contained in the 
continuous action, the discrete symmetries only come into play when F3 is rational. 

An element (/ci, A;2, k^) G K acts on Aq as (ai + /ci, 02 + A;2, ^3 + ^3), and so it acts in 
coordinates as 

{Z + {ki-k2),Y + ki). (3.10) 
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But note that {ki — ^2) = ^3{k2 — ks), and thus the action of K on 7? just adds Fa times 
an integer to the coordinate Z. When is irrational, ki = k2 = ks, and K does not act on 
the Z component. When Fs = p/q is rational (in lowest terms), ki — k2 is actually p times 
an integer because the ki are integers. Thus the induced symmetry on R is translation by 
p. Algebraically speaking, the map {ki,k2,k^) 1— >■ {ki — k2,ki) is an isomorphism from K 
to pZ © Z. In this new form K acts on R with just its first factor. Thus the action of K 
on R induces a mod p symmetry on the plane R. The Hamiltonian on R thus has a mod 
p symmetry and so the final factored system has a cylinder of width p as its phase space. 
Figure 3.2 shows the orbits of the reduced system for a vortex array with Fi = 1,F2 = 1/2 
and F3 = —3/2. The left and right edges of the box can evidentially be identified to get a 
width 3 cylinder. 

In the case of lattices, there is a discrete action on both the real and imaginary parts of 
Za- Thus if F3 = p/g, the factored system may be considered a torus with width and length 
p, i.e. a Hamiltonian on M^/pZ^. 

§3.5 Integral manifolds and dynamics. As noted in §2.2, in the case N — 3 the vortex 
systems arc completely integrable as Hamiltonian systems. Orbits of the phase flow are 
constrained to simultaneous level sets of H, Q and P. If the phase space was compact, these 
integral manifolds would be generically tori. The systems here do not have compact phase 
spaces, but the reduction and subsequent reconstruction make the structure of the integral 
manifolds clear. 

In the reduced space R the generic, bounded level set of Kx is a topological circle 
(dynamically a periodic orbit). If we fix one such circle C G R, then an integral manifold A 
of the full system corresponding to fixing values of H, Q and P is the preimage of C under 
the projection Aq ^ R. Thus the generic A is homeomorphic to C x with the factor 
coordinatized by Y. The fiow on A factors as in (3.8). The simplest way to understand this 
fiow is to examine the return map of orbits to the M? factor as one travels around C. By 
(3.9) and (3.10) this return map is translation by b{ZQ) for any Zq E C. 

In the array case when F3 is rational we have used the action of the discrete symmetry 
in the first coordinate of (3.10) to factor R. One also sees from (3.10) that the action of 
the discrete symmetry on the Y-component is addition of an integer. This action can be 
used to factor the second component of A to obtain an integral manifold homeomorphic to 
C X (S^ X R) with the return map to the cylinder factor now addition of 6(Zo) (reduced 
mod 1 in the first coordinate). The lattice case is similar, but the discrete symmetries yield 
integral manifolds C x (S^ x S^). Depending on the value of 6(^0), the fiow on these tori 
can have 1, 2 or 3 independent frequencies. 

§3.6 Irrational circulation ratios and quasi-periodicity. We now use the continuous 
and discrete symmetries from a slightly different point of view with the goal of illuminating 
the case of irrational circulation ratios and the resulting quasi-periodicity in the reduced 
Hamiltonian system. We factor by the action of the discrete and then the continuous sym- 
metry. The level sets of J descend to invariant subsets, but they are no longer level sets of 
a real-valued function. The focus here is on the lattice case and we comment on the array 
case at the end. The nomenclature n-torus refers to the topological space T"" := {S^)'"', i.e. 
the n-dimensional manifold that is periodic in all n directions. 

Recall that the phase space of the Hamiltonian system describing the vortex motions 
is — T^, which is denoted here as y. The discrete action of Z^ on y factors the phase 
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space to {S^ x S^)^ — T' which wc give coordinates (i^i, ^2, i^s), where T' is the factored 
coUision set. (Note that each x S*^ is a Lie Group; one does arithmetic by reducing mod 
1 in both factors.) The continuous symmetry of acts by simultaneous addition on all 
2-tori, so T e acts as {ui + r, i/2 + r, us + t). Identifying points under this action yields 
a quotient space y — {S^ x S^)'^ — T with coordinates {ui — V2, 1^2 — v^) '■— {U, V) and 
t = {[/ = 0} U {y = 0} U {[/ = V}. The Hamiltonian (2.4) on y descends to 

H{u, V) = -^(rir20(c/) + V2Vi^{v) + rir30(c/ + v)). (3.11) 

If is an integer, the integral J may be expressed in these coordinates as J{U, V) = U — TojV 
(recall that Fi = 1) with J having values in a 2-torus with width F3 in both directions. If 
Fs is irrational, no such adaptation may be made. But in any case, since a level set Aq of 
J is invariant under the Hamiltonian flow, the projection Aq to y is invariant under the 
Hamiltonian flow induced by H. 

We now treat j) as a subset of the 4-torus, T"^. If F3 is irrational, the level set ^0 is 
an immersed two-plane that winds densely in the four torus. The Hamiltonian H on 
restricted to is the same as the reduced Hamiltonian on R given in (3.6). Thus this 
Hamiltonian induces a quasiperiodic system on R in the sense of [Al]. The function H has 
singularities when U = 0, V = OorU = —V, each of which represents a plane in T^. The 
intersection of these planes with the densely wrapped plane Aq yield the poles of the reduced 
Hamiltonian K£. Thus the collection of poles are a quasi-crystal, at least according to one 
deflnition of that term ([J], [Al]). 

Since (the complexiflcation of) H has simple poles, the poles of the restricted Hamiltonian 
on any Aq will also be simple. Thus it is plausible that the flow on R (or Aq) can be treated 
as advection in the presence of fixed vortices with circulations given by the residues of the 
poles. The precise sense in which this is true is given in [SA] and [AS2] . Further, it is shown 
in [AS2], that this collection of fixed vortices is dynamically fixed as well, i.e. if these vortices 
are allowed to interact with each other, they will still be stationary. In other language, the 
collection of poles of the reduced Hamiltonian on R forms an equilibrium configuration of 
the planar A^-vortex problem. 

This viewpoint can also be used to illuminate what happens as the various parameters 
in the system are varied. For example, once the F^ are fixed, changing Jq just translates 
Aq and does not alter its "slope". Thus there will be no qualitative changes in the pattern 
of the poles of Aq except for the exceptional values of Jq which cause it to hit the triple 
intersection points U — V — (the integer lattice). On the other hand, if the F^'s are varied, 
the background fiow changes as does the slope of the Aq. However, the singularity planes 
remain the same, so it is easy to track the singularities on Aq. In particular, if there is a 

sequence of rational Fg"'-' = ^ with ^ ^ ^3 with F3 irrational, then the periodic systems 
will converge to the quasiperiodic one after the appropriate rescalings. 

A similar analysis can be applied to the array case, but now if F3 is irrational the quasi- 
periodicity of the reduced Hamiltonian Kj\^ will only be in the real direction. 

§4 Regimes and braids. 

Braids are the standard mathematical tool for describing the topology of periodic motions 
of collections of points in the plane. In this section braids are used to describe the motions 
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of point vortices. While a braid may be used to describe the periodic motion of any number 
of vortices, the analysis is particularly simple in the case of three vortices. As we have seen 
in §3 when the total circulation is zero, the generic motion in the reduced plane is periodic. 
Thus the corresponding motion of the vortices is periodic after adjusting for the dynamic 
phase, and so may be described by a braid. 

§4.1 Braids and braid types. This subsection gives a brief introduction to braids from a 
point of view useful for the applications that follow. For broader perspective and additional 
information sec [Bml] or [BL]. Braids are used to describe the motions of the bodies in the 
planar A^-body problem by Montgomery [Mt]. 

A physical or geometric braid is a collection of noncrossing paths in that start at some 
finite collection of points E on the plane where the third coordinate is zero and end at the 
same set of points (prhaps permuted) on the plane where the third coordinate is one. More 
formally, a physical braid on n strands is a collection of maps b = {bi, ... , with each 
bi : [0, 1] and (1) each bi has the form bi{t) = (a.j(t),t) with each continuous, (2) for 

all t G [0, 1], bi{t) ^ bj{t) for i ^ j., and (3) the set of initial points is identical to the set 
of final points E — {ai(0), 02(0), . . . , 0^(0)} = 02(1), ■ ■ ■ , The initial and 

final points of the braid are collectively called the endpoints. The collection of trajectories 
{ai{t)} is called the plane motion of the braid. 

A mathematical braid on n strands is an element of the braid group Bn- This group is 
defined as possessing the n generators ai,a2, ... , cr„ with inverses denoted ai,a2, ... , a"n 
and the relations aiaj — UjUi for \i — j\ > 1 and (JiUi^iai — ai+iaiai+i for all i ([Bml]). 
A braid word refers to a sequence of "letters" , with each letter being one of the a-i or their 
inverses. The inverse of a generator is indicated by an overbar, CTj. Two braid words are 
said to be equivalent if they represent the same element in the braid group, i.e. one can be 
transformed into the other using the relations in the group. The identity element in Bn is 
denoted e. 

The assignment of a mathematical braid to a physical one requires a plane onto which 
projections arc made. The convention here is to let this plane in he y = k for some 
large negative k. If we treat the xy plane as the complex numbers, projection of b onto 
the chosen plane yields a family of curves (Re(aj(t)), t). An ij-crossing is a point where 
Rc{a.j{t)) = Re{aj{t)) for some t and i ^ j. The braid word will record which strand is in 
front at each crossing or, more precisely, which of the 6j(t) is closest to the projection plane. 

The braid word corresponding to b is read off from the picture of the projection. Assume 
for the moment that all crossings are transverse and take place at distinct times. Starting 
from the bottom, examine the first crossing. If the i*^ strand from the left crosses behind the 
(i + 1)**, write down the letter (jj. If it crosses in front, write ct^. Now continue upward in the 
projection checking each crossing and writing a braid letter. Note that the i^^^ strand at each 
step refers to the i^^ strand from the left at that level. The physical strand that is i^^'' at one 
level may become (i + 1)*** or {i — 1)*** at the next level. It is also worth noting that contrary 
to the conventions of [Bml] , the convention here is that letters further to the right in a braid 
word encode crossings that are higher up the braid. We adopt this somewhat nonstandard 
convention because our physical braids arise from time-parametcriscd trajectories in the 
plane, and it is more natural to visualise these in with time going upwards. 

It is natural to ask how much the mathematical braid tells us about the physical braid. 
In particular, which physical braids get assigned the same mathematical one? Say that two 
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physical braids b and b' with the same endpoints are equivalent if one can be obtained from 
the other by a deformation that fixes the endpoints and does not cut the strands, i.e. there 
is a continuous family of physical braids b^ for s G [0, 1] with b*^ = b and b^ = b'. The 
relations in the braid group are chosen precisely so that two physical braids are equivalent 
exactly when they are assigned equivalent braid words, i.e. the same mathematical braid. 
A theorem of Artin's ([Bml]) says that any equivalence between physical braids can be 
described in terms of just the two kinds of deformations described by the relations in the 
braid group. 

Because equivalent physical braids are assigned equivalent braid words we can eliminate 
the assumption we had to make in order to assign a mathematical braid to a physical one. 
If b does not have transverse and time-distinct crossings, deform it to another braid b' that 
does and compute a mathematical braid for b'. Since any other "good" deformation will be 
assigned a braid word equivalent to that of b', this braid word may be used unambiguously 
for b. 

Thus far we have restricted attention to equivalence of physical braids with the same 
endpoints. Since the braids here arise from vortex motions with a variety of initial positions 
we need to extend the notion of equivalence. Informally two braids with perhaps different 
endpoints are equivalent if one can be transformed into the other via a plane transformation 
applied on each horizontal plane followed by a deformation with fixed endpoints. More 
formally, two physical n-braids b and b' are equivalent if there is a homeomorphism /i : — > 
with h{ai{0)) = a[{0) for all i and the physical braid {{h{ai{t)),t)} is equivalent with fixed 
endpoints to b'. A simple argument shows that physical braids that are equivalent in this 
sense are assigned conjugate elements of the braid group, i.e. the words satisfy w' = gwg~^, 
with g an element of the braid group that represents the homeomorphism h. Note that 
changing the projection plane may be accomplished by a rotation in the plane, and so the 
mathematical braid obtained using the new projection plane will be conjugate to the original. 

Since the goal here is to use braids to analyse the topology of vortex motions, topolog- 
ically equivalent physical braids and physical braids viewed by different observers must be 
assigned the same mathematical object. The remarks of the last paragraph make it clear 
that this object cannot be just a mathematical braid, but rather must be a conjugacy class 
in the braid group, i.e. the collection of all elements conjugate to a given one (and thus 
conjugate to each other). These conjugacy classes were called braid types in a related context 
and that terminology is adopted here (cf. [Bd]). For simplicity of exposition in the sequel 
we will often speak informally of the braid associated to a physical braid, but a more careful 
terminology would be braid type. 

§4.2 Three vortices in the plane. In assigning braids to the motion of three point vortices 
in the infinite plane we continue to focus on the case of zero net circulation, Tq = 0. As 
in §3, we fix a value of the integral Jq = Qo + iPo, perform the reduction, and obtain the 
Hamiltonian (3.6) on a copy of the plane R with complex coordinate Z. 

Pick an initial condition Zq in R which is contained in a periodic orbit Z{t) := Z{t] Zq) 
with period P. Using (3.4), the motion Z{t) determines the evolution of the differences in 
the positions of the three vortices. Thus Z{t) can be used to compute the motion in the 
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frame of one of the vortices. In the 2;2-frame this motion is described by 
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° (4.1) 

Jo - riZ(0 

Thus all the ZaS are periodic with period P. Note that Za ^ zp lor a ^ (3 (i.e. the absence 
of collisions between vortices) is equivalent to Zq, ^ Zp for a ^ p. Again using (3.4), this 
happens as long as Z{t) avoids the points pi := Jo/^i,P2 '■— —Jq/^2, and ps :— 0. The 
points Pa are the poles of the Hamiltonian K-p from (3.6). 

The vortex motions generate a physical braid by treating the vertical direction as time 
and defining 

Za^{Za{t),t), (4.2) 

for a = 1, 2, 3. Each path Za then connects a point on the plane t = to the same point on 
the plane t = P, and further, since Z{t) does not pass through any of the poles, the paths do 
not intersect. Thus the paths yield a physical braid on three strands, which may be assigned 
an element of as in the previous subsection. 

This assignment is most easily accomplished by taking a more topological point of view. 
If we define R" = R— {j5i,P2iP3}i then the periodic orbit Z{t) in R represents a closed curve 
(or loop) in R° . The procedure used to construct a physical braid from Z(t) can equally 
well be used to assign one to any loop in R° . If 7 is a loop in R°, i.e. 7 : [0, 1] — > R° with 
7(0) = 7(1), let Zi, Z2, Z3 represent the three paths given by (4.1) using Z{t) — ^{t). The 
physical braid corresponding to 7 is then generated by (4.2). 

The computation of the mathematical braid corresponding to this physical braid just 
requires knowledge of the crossing of its strands. Since we are viewing from the negative 
imaginary axis, the positions of the Za from left to right are determined by their real parts. 
Strands cross exactly when this order changes. This can only happen when 'Re{Za{t)) — 
Re(ZQ,/(t)), or equivalently, when Re(Z(t)) = Re(p^) where (3 is the index different from a 
and a' . Thus the vertical lines through the poles divide R° into vertical strips in which the 
order of the Za from left to right is constant. These vertical lines will be called crossing lines. 

To specify the generator corresponding to the crossing of one of these lines it is also 
necessary to know which strand is in front. This information is determined by the sign of 
lm{Z{tc)) — Im(p^) when Z{tc) lies on the crossing line Re(Z) = Re(p^). This sign can only 
change on a crossing line when Im(Z(tc)) — Ini(p^) = 0, i.e. at the pole. Thus all the crossings 
on the same side of the pole correspond to the same strand being in front. We will call a 
component of a crossing line minus its pole a generator arc. The final piece of information 
needed to specify the generator describing a crossing is the direction that 7 traverses the 
generator arc. In summary then, any two crossings of a generator arc in the same direction 
contribute the same generator to the braid. Thus to compute braids of loops or actual vortex 
motions it suffices to compute the generators corresponding to each generator arc. A sample 
computation is given in the next subsection. 

There are two issues that need to be clarified before moving on to the example. First, 
unlike a mathematical loop there is no distinguished starting point for a periodic motion of 
the vortices. However, changing the initial point on Z{t) will only cyclically permute the 
generators in the braid word representing the motion. This corresponds to conjugating the 
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word in the braid group and this ambiguity has ahcady been resolved using the braid type. 
The second issue is that the usual generators of the braid group only keep track of which 
adjacent strands are crossing and not the numbering of the stands. Thus in assigning a 
braid generator to a crossing we have lost track of which vortex has crossed which; our braid 
word is only encoding the topological type of the interaction, and not whether, say, vortices 
2 and 3 are the pair that are circling each other. However this information can be recovered 
easily by knowing the relative positions of the vortices at the starting point of the braid and 
following the strands. 

§4.3 A planar example. Figure 4.1 shows the crossing lines for the planar case of Figure 
3.1. The pair of numbers at the top of the crossing line indicate which pair of vortices are 
crossing. Each generator arc has a number j and an arrow; crossing the arc in the direction 
of the arrow contributes the positive generator aj to the braid word. Crossing in the opposite 
direction contributes the inverse of this generator aj. The braid word corresponding to a loop 
7(t) can be computed by writing in order the generators corresponding to 7(t)'s sequential 
crossings of generator arcs. 




Figure 4.1: The crossing lines and generator arcs for the system 
shown in Figure 3.1. 

As an example, let us compute the braid word corresponding to periodic orbit Z'{t) 
indicated in the region labelled D in Figure 4.1. Begin tracking the motion at the starting 
point labelled Z'{tQ) in the upper left hand corner. At this point the vortices from left to 
right are 1,3,2. The first generator arc crossed is above the pole p2, it contributes a ai to 
the braid word. Moving across horizontally, the next generator arc contributes a (T2, and so 
on. When we return to the initial point the entire braid word aia2(Ji(Jio'2(Ji has been read 
off. This braid is shown in Figure 4.2c, with the strands labelled below. By examining the 
braid we see that this motion corresponds to the vortices rotating clockwise around each 
other. Examining the braid one can see that the motion is topologically the same as all the 
vortices rotating once around a circle clockwise. This becomes especially clear after using 
the braid relation and rewriting the braid as ((Ti(72)^. 

Figure 4.2a shows the planar trajectories of three vortices corresponding to the periodic 
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Figure 4.2: Vortex motions from regime D of Figure 3.1. (a) Tra- 
jectories of the three vortices in the plane, (b) Physical braid of this 
motion in the frame of the second vortex, (c) Mathematical braid of 
the motion. 

orbit Z'{t). Recall from §3.2 that there are many vortex motions corresponding to a given 
periodic orbit in R, but all these motions differ by a uniform translation, so it suffices to 
pick an initial position of one of the vortices, say 2:2(0) = 0. Figure 4.2b shows the physical 
braid corresponding to the vortex motion in the frame of the second vortex as given by (4.1) 
and (4.2). 

The other regions A, B, and C each contain a pole. The braid describing each of these 
regimes is the square of a generator and corresponds to a pair of vortices rotating about 
each other once per period. There is no linking with the third vortex. The pair of vortices 
that are involved in the rotation depends on the pole in question, and the direction of the 
rotation (i.e. whether the braid is af or af) depends on the signs of the circulations of the 
interacting vortices. 

Recall that a separatrix is an orbit connecting two saddle points. Define a regime as a 
connected component of R— {poles and separatriccs}. It is clear from Figure 4.1 that every 
periodic orbit in the same regime yields the same sequence of crossings of generator arcs, 
and thus the same braid word. Therefore all the motions within a regime have the same 
topological type of motion. More generally, any pair of loops that can be continuously de- 
formed into each other (i.e. are homotopic) in R° will be assigned the same braid; homotopic 
loops may cross different generator arcs, but any extra crossings will consist of a collection 
of crossings and then reverse crossings of the same generator arcs. From another point of 
view, deforming the loop corresponds to deforming the physical braid. The corresponding 
mathematical braid will only change when a pair of strands go though each other. This 
can only happen when the loop passes through a pole, which is not allowed as all our loop 
deformations are in the complement of the poles. 

This situation is somewhat analogous to the residue theorem in which a deformation 
of a closed path in the complement of the poles does not change the value of the integral. 
However, the situation here is more restrictive; homologous loops yield the same integral but 
only homotopic loops give the same braid. In algebraic language, after fixing a basepoint for 
loops, the process described here gives a homomorphism from the fundamental group of the 
punctured plane, 7ri(i?°), to the braid group on three strands, S3. 

Although all periodic orbits in the same regime have the same braid, it is important to 
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note that the period P and dynamic phase b of the periodic orbits in a single regime can 
vary greatly. The period will go to zero near a pole and approach infinity as orbits near a 
separatrix. Thus the time scale of the actual vortex motions in the same regimes can be very 
different. In addition, variations in the dynamic phase can make a significant difference in 
the motion as observed in the lab frame. 

§4.4 Three-vortex arrays. This subsection develops the tools needed to assign a braid 
to the periodic motion of three vortices in a singly-periodic array. We continue to restrict 
to the case of zero net circulation X! Ta = 0, fix a value of the integral Jo = Qo + 'iPo, and 
perform the reduction as in §3. The result is a Hamiltonian system given by (3.6) on a copy 
of the plane R with complex coordinate Z. If Fa is the rational p/q (in lowest terms), then 
the system on it! has a mod p symmetry in the real part. 

By treating each vortex position 2q dbS db sill gly-periodic variable, the vortex motion can 
be viewed as taking place on a cylinder x R (the circle here has perimeter one since 
we continue to restrict to the case L — 1). Equivalently, in the language of §3, we factor out 
all the motions by the discrete symmetry. Again we eliminate the dynamic phase by passing 
to the frame of the second vortex, thus making the motion periodic. Accordingly, let 
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in which properly speaking the subtraction is done in the Lie group x R. The triple 
(ci(t), C2(t), cs{t)) is called the cylinder motion of the vortices. 

A braid describes the motion of a set of points in the plane. To translate the cylinder 
motion to the plane, recall that the cylinder is topologically equivalent to the punctured 
complex plane C — {0}. This topological equivalence is realized by the conformal map 
T[z) = exp{27riz). Under this conformal map a path around the cylinder transforms to a 
path around the origin in the complex plane. To capture this type of motion in the braid 
description, we add the constant path at as the last coordinate of the motion in the plane. 
Thus the plane motion of the vortices is given by 

si(i) := T{ci{t)) = exp{27riZ{t)) 
S2{t) :=r(c2(t)) = l 

m T( (9 - 'Ai-riZlf) , (4.4) 

ssit) ■= T[c2{t)) = exp{2m ) 

i 3 

S4it) = 0. 

Using the same technique as in the previous subsection we associate a braid (now on 4 
strands) to the plane motion of a three-vortex array. Again we visualise the plane motion 

in three dimensions using (4.2) with time going upward and view the resulting strands from 
the negative imaginary axis. The crossing curves in R describe the positions of an orbit 
Z{t) at which Sa and 5^3 cross; this happens when their real parts are equal. To express 
this in equation form, define four functions on R by ipi{Z) — T{Z),il)2{Z) — l,ips{Z) — 
r((Jo - riZ)/r3), and MZ) = O. Then s«(i) = MZit)) = T{cait)), and the {a, (3) 
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crossing lines are defined by the equation Rc{ilja{Z)) = Re{ipi3{Z)). In contrast to the planar 
case, these equations no longer yield just straight lines and thus the change in terminology 
to crossing curves. 

The strand which is in front at a crossing is determined by the sign of lm{sa(tc)) — 
lm{sp{tc)) when 7(tc) is on the corresponding crossing curve. This sign changes on the 
crossing curve exactly when Sa = sp, which corresponds to the poles of the Hamiltonian 
(3.6). Note that it is no longer the case that every crossing curve contains a pole. This 
means that crossing anywhere on the curve yields the same generator of the braid word, i.e. 
the entire curve is a single generator arc. 

§4.5 An cirray example. We now perform the calculations described in §4.4 on the array 
case shown in Figure 3.2. See [SA] for a details on this and other array examples. The 
crossing curves corresponding to si and S3 are defined by Re{ipi{Z)) — Re{ips{Z)), which, 
writing Z = x + iy, is 

2 1 2 1 

exp(-27ry) cos(27rx) = exp(-27r(-?/ - -)) cos(27r(-a; - — )). 

Both cosine terms can vanish yielding the two vertical lines x — 5/4 and x — 11/4. Otherwise 
we may solve for 

3 , cos(27ra;) , 

which is defined only on the (mod 3) intervals (-3/4, 1/4), (1/2, 3/4), and (7/4,2). Thus 
there are five distinct (1, 3)-crossing curves. The computations for other pairs is similar and 
the various crossing curves are shown in Figure 4.3a. The numbers separated by a hyphen 
indicate which pair of s^s are crossing. Note that since the vertical lines x = 5/4 and 
X = 11/4 correspond to = Re(si) = Re(s3), and we are considering as the position of S4, 
these vertical lines correspond to a triple crossing. 

The poles are the points 0, 1 and 2 corresponding to si = S2, the points | + |^ and 
+ |i corresponding to S2 = S3, and the point —j — ji corresponding to si = S3. Note that 
many crossing curves do not contain a pole. Figure 4.3b shows the results of computing the 
generator arcs, with labels using the same conventions as in Figure 4.1. 

It is important to note that the details of Figure 4.3 are not intrinsic to the vortex 
motions. Changing the projection plane used to compute the braids would change all the 
crossing curves. However, the braid which is computed for a loop is intrinsic up to conjugacy. 
In algebraic language, after fixing a basepoint, the equations (4.4) (using Z(t) to represent 
a general loop) may be used to define a homomorphism of 7ri(i? — {poles}) B4. This 
homomorphism is not intrinsic, but it is intrinsic up to conjugacy in i?4. 

Now define a regime as in the planar case, i.e. as a connected component of i? — 
{poles and separatrices}. Each periodic orbit in a regime gives rise to the same mathe- 
matical braid; these are listed in the second column of Table 4.1. Note that the relations in 
the braid group have been used to simplify many of the braid words. The braids for regimes 
I and IX through XIII were computed using closed trajectories with initial position near the 
point labelled r in Figure 4.3b. The initial position used for regime II was near the point 
labelled q. 

As a sample we compute the mathematical braids corresponding to regimes I and XIII. 
The closed orbits in regime I travel across the top of the reduced plane from left to right. To 
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Figure 4.3: (a) Crossing curves and (b) generator arcs for the sys- 
tem in Figure 3.2 (the picture in (b) has been distorted for increased 
clarity). The points labeUed q and r correspond to initial points of 
trajectories as explained in the text. 



compute the braid we chose the point labeled r in Figure 4.3b as the initial point. Moving 
to the right from this point we first traverse a vertical line on which there is an arrow 
labeled 1. This indicates that the motion of the trajectory has resulted in a crossing of 
the two left most strands in the physical braid representing the vortices. Further, since the 
trajectory's motion is opposite to that of the arrow, the strands are crossing with the left 
one in front of the right one. Thus the traversing of the first crossing curve contributes 
the letter a\ to the braid word of the trajectory. Continuing to the right the trajectory 
traverses a crossing curve on which there us an arrow labeled 2 and the arrow points in 
the direction opposite to the traverse. Thus this motion contributes a (72 to the braid word 
of the trajectory. Continuing across the top of Figure 4.3b, we add braid letters cxi, then 
a"2, etc. After encountering the right edge of the phase plane the identification with the 
left edge allows a return to the initial point. The entire braid word for this regime is thus 
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Regime 


Braid 


TN type 


Expansion 
Constant 


I 


(aiCT2)^CT?(a-25-i)3 


reducible, all f.o. 




II 


0'3f2(o'lCr2(T3)^cr2Cr3Crif72a-i((TitT2Cr3)^CT2<TiC73CriCr2(Ti(T3 


reducible, all f.o. 




III 


4 


reducible, all f.o. 




IV 


4 


reducible, all f.o. 




V 


4 


reducible, all f.o. 




VI 


-4 


reducible, all f.o. 




VII 


4 


reducible, all f.o. 




VIII 


4 


reducible, all f.o. 




IX 


1 2 

44 


reducible, all f.o. 




X 


(<73CT2)^ 


reducible, all f.o. 




XI 


1 \2 2 /- - \3- 
lO'3(T2j 0'30'lO'20'iO'20'lO'20'30'l'''3lO"lO"2j cri 


pA 


13.93 


XII 


{oz02)^ozo\a2o\{a20zo\az)^(p\a2fo\ 


reducible, one pA 


13.93 


XIII 


(aiCT2)^CT|(CT2a-i)^ 


pA 


9.90 



Table 4.1: The mathematical braids and Thurston-Nielsen type of 
the regimes for the system shown in Figure 3.2. The expansion con- 
stants for pseudoAnosov (pA) components are given. In the finite 
order (f.o.) case these constants are 1. This means that there is no 
intrinsic topological expansion in that case. 



^1^2^1^2^1^2^1^1^2^1^2^1^2^1 = (^1^2)^^l(^2^l)^- 

The closed trajectories for regime XIII have a motion similar to those of regime I with 
the crucial topological difference that they pass below the pole p3 which is located at 1 + Oi. 
This implies that in the complement of the poles, the closed loops corresponding to regimes 
I and XIII are not homotopic. Beginning a closed trajectory representing regime XIII near 
the point r, the braid word starts with G\a20\a2 like that of regime I, but since the loops 
for regime Xlll pass below the point p3, their braid words have a a^oz next, and then 
continue with the same letters as regime I. The entire braid word for regime XII is thus 
0\a20\a20z(yz(y\020\G\G20\a20\a20 Note that the braid relation a\a2Cf\ = a2^i&2 used 
on the triple of letters right after the a^^a-i allows us to reduce the word to {oia2)'^o-1{a20i)^ 
as given in Table 4.1. 

Wc now briefly describe some of the motions. In these descriptions we will use the 
word "vortex" to refer to the trajectories of the planar motions, si,S2,S3, and 54. Properly 
speaking si, S2, and S3 are the transformed positions of the vortices, and S4 is not a vortex 
at all, but rather it keeps track of the position of the origin. 

Before interpreting the braid, it is useful to recall the various transformations that have 
been used. The motion in the singly periodic plane is really of infinite families of vortices. 
These famihes are treated as three individual vortices on a cylinder and then this motion is 
transformed to the punctured plane. Thus a strand in the braid rotating about the strand 
representing vortex 4 (= the origin) describes motion around the cylinder, which in turn, is 
a horizontal translation by one period of the corresponding family of vortices in the singly- 
periodic plane. In addition, the braid is obtained using the frame of the second vortex, so in 
particular vortex 2 will never rotate around the origin, and thus the strands corresponding 
to vortex 2 and 4 will always be parallel. 
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The regimes with the simplest braids are those that contain a pole, III through VIII. The 
vortex motions are topologically the same as regions A, B, and C in the plane case: a pair 
of vortices are circling each other and the other vortex is uninvolved. There is no rotation 
around the origin, thus there is no net rotation around the cylinder. 

The vortex motions corresponding to regime X also have no net motion around the 
cylinder (see figure 6.1a). This is indicated by the fact that none of the other strands link 
with the far left one which represents the origin of C. Examining the right hand sub-braid 
on 3 strands, one sees that every vortex rotates about every other one once in the clockwise 
direction. The net motion is the same as that in the planar region D and it amounts to 
one full rotation of the vortices as a group. The motions in regime IX are similar, but the 
sub-braid is slightly different; vortices 2 and 3 rotate about each other clockwise, and then 
vortices 1 and 3 do likewise, but vortices 1 and 2 do not link. 

The regimes with the next simplest braids are I and II. The mathematical braid for I is 
shown in Figure 4.4d. Vortex 1 rotates three times around the origin counter clockwise and 
vortex 3 rotates twice, also counter clockwise. This rotation takes place "inside" the fixed 
vortex 2. In the singly periodic plane, this corresponds to the vortex 1 family translating 
three strips to the right and the vortex 3 family two strips; this translation takes place 
below the position of vortex 2. This motion is shown in Figure 4.4a with one vortex from 
each singly periodic family pictured. Figure 4.4b shows the motion in the frame of the 
second vortex after it has been transformed to the plane. The physical braid of this motion 
constructed using (4.4) and (4.2) is shown in Figure 4.4c. The motion in regime II is similar, 
but now the rotation is clockwise, and goes around vortex 2, which simply means that in the 
singly-periodic plane the translations of vortex 1 and 3 are above the position of vortex 2. 

The motion in regimes XI, XII, and XIII is sufficiently complicated that the braid itself 
is the best description of the motion. One feature of interest in regime XII (see Figure 6.1b) 
is that the vortex pair 1 and 3 have the same motion with respect to the other two strands 
while they rotate about each other three times. This kind of hierarchy of motion is described 
precisely by the notion of reducibility introduced in §5. 

§4.6 Three-vortex lattices. The method of analysis of this section can also be apphed to 
vortex lattices. When Fs = p/q is rational, the reduced Hamiltonian system (3.6) can be 
viewed as taking place on a torus with width p in both directions. The generic orbit of this 
system will be a periodic orbit Z{t). This periodic orbit can be used to generate a motion 
of the three vortices on the torus which will be periodic in the frame of one of the vortices. 
However, in contrast to the array case, this cannot be transformed into a motion on the 
plane, and so one must record the motion using the braid group of the torus. This group 
not only has generators corresponding to the switching of two points as in the usual braid 
group, but it also has generators describing motions that circulate around the meridian and 
longitude of the torus (see [Bm2]). Thus the braids are more difficult to visualise. However 
the basic features of the analysis in the array and plane case go through with no difficulty. 
In particular, each regime is assigned just one braid, and so the braid captures the sense 
in which all the vortex motions in one regime are topologically the same and distinct from 
those of other regimes. 

§5 Isotopy classes and the Thurston-Nielsen theory. 

This section provides an introduction to isotopies and the Thurston-Nielsen theory in 
preparation for the study of the advection homeomorphisms generated by vortex motions. 
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Figure 4.4: Vortex motion from regime II of Figure 3.2. (a) Tra- 
jectories of the three vortices in the singly-periodic plane (only one 
representative of each family is shown) . (b) The motion in the frame 
of the second vortex transformed into the plane, (c) Physical braid of 
the transformed motion, (d) Mathematical braid of the transformed 
motion. 

Roughly speaking, two homeomorphisms are isotopic if one is a continuous deformation of 
the other. The braid of a vortex motion will be seen to determine the isotopy class of the 
corresponding advection homcomorphism. The braid is then used as input into the Thurston- 
Nielsen theory which describes dynamical data that is shared by isotopic homeomorphisms, 
and is thus present in the advection homcomorphism. 

§5.1 Homeomorphisms, isotopy, and braids. A one-to-one and onto map f : X ^ X is 
called a homeomorphism if / and its inverse, / \ are both continuous. If in addition, / and 

are both differentiable, then / is called a difFeomorphism. Two homeomorphisms / and 
g are conjugate if there is a third homeomorphism h with / = hgh~^. The homcomorphism 
h represents the change of coordinates from / to g. If £^ is a finite set and f{E) = E, then 
/ is called a rel E difFeomorphism or homeomorphism. 

Most of the homeomorphisms occuring in this paper are defined using the advection of 
a passive particle in a time varying velocity field. Given a vector field on the plane {u, v) , 
the advection equations are 



dx dy . , 

-^ = u{x,y,t) -^ = v{x,y,t). 



(5.1) 



Letting z = {x,y), the solution with position zq at time to is denoted z{t; -ZQi^o)- By con- 
vention the running time t satisfies z(tQ; zqAq) = zq, not z{0; zo,tQ) = zq (unless, of course, 
to = 0). Fixing a time t, the time t-flow map or the time t-advection homeomorphism is 
defined as ht{zo) := z{t;zo,0), i.e. we advect each particle for time t, and ht maps each 
initial position to its final position. The collection of ht for all t describes all the solutions 



23 



and is called a fluid motion (cf. §3.1 on the use of "flow"). Note that given a fluid motion 
the underlying velocity fleld is obtained by difl^erentiation. 

The usual deflnition of isotopy rel a flnite set E is the following. Given two rel E 
homeomorphisms /j : ^ M^, say that /o is isotopic to fi rel E if they can be connected 
by a continuous family of rel E homeomorphisms ft, t e [0, 1]. Since £^ is a flnite set this 
implies that ft{k) = fo{k) for all t and k E E. The collection of homeomorphisms isotopic 
to / rel E is called its rel E isotopy class. The flnite set E in the sequel will be the initial 
positions of the collection of point vortices. These points will be flxed under the advcction 
homeomorphism considered, but in general, the points of E could be permuted by a rel E 
homeomorphism. It is sometimes convenient to remove the set E from the plane and consider 
a rel E homeomorphism as a homeomorphism of the punctured plane — E. Note that two 
rel E homeomorphisms are isotopic rel E exactly when they are isotopic as homeomorphisms 
of m2 _ E. 

The exposition of isotopy in fluid mechanical terms is most clear in the special case of 
isotopy to the identity rel E. For simplicity we restrict to the case of diffeomorphisms as 
is natural for fluid mechanical applications. A rel E diffcomorphism / : — ^ will be 
isotopic to the identity rel E exactly when / is the time 1 advection homeomorphism of a 
velocity fleld which satisfles u{k, t) — v{k, t) — Q for all /c e £■ and t e [0, 1]. This means that 
all the points of E are flxed throughout the advcction. Thus we may think of / as arising 
from advection due to body forces on the fluid and the points in E as flxed point-like regions 
that are not permeable by the fluid. If the points of E move under a given fluid motion, the 
time-one flow map can still be isotopic to the identity if there is another fluid motion that 
also gives rise to /, but in this new motion the points of E are flxed for the entire motion. 

The braid associated with the motion of points in E can be used to determine whether 
a given advection diffcomorphism is isotopic to the identity. Given a fluid motion ht^ the 
motion of the points in E is given by aiit) := ht{ki) for each ki G E. Now a family 
of diffeomorphisms V't : ^ for t e [0, 1] with ipt{ai{t)) = ki for aU ki e E and 
V'O = "^1 = id induces a level preserving diffcomorphism ^ : x [0, 1] — >■ x [0, 1] dcflned 
by "^{Zjt) = {ipt{z),t). The diffcomorphism ^ will take b to the trivial physical braid which 
consists of just vertical strands at the points of E. If this is the case say that b is weakly 
equivalent to the trivial braid htriv 

The relationship of braids to isotopy classes is this: given a fluid motion ht with advection 
homeomorphism f = hi and a flnite set E with f{E) = E, then / is isotopic to the identity 
rel E if and only if the physical braid b describing the motion of E is weakly equivalent to 
hfriv Here is the proof. If b is weakly equivalent to htriv by the family ipt, then h'^ = ipt^ ht 
is a fluid motion that flxes E for alH e [0, 1] and h[ — hi — f, and so / is isotopic to the 
identity. Conversely, if h't is a fluid motion that flxes E for all t G [0, 1] and h'l — hi — f, 
then ipt = h[o h^ is a family that will make b weakly equivalent to \)tfiv- 

The next step is to see how the braid word describing the motion of E indicates when 
the advection diffeomorphism is isotopic to the identity rel E. This is certainly the case 
when the mathematical braid is the identity element in the braid group, but this is not the 
only case. Consider the plane motion given in complex coordinates by ht{z) = \z\ exp(27r?'t). 
For a given n, let E be the set of n*'' roots of unity, {exp 27ri j'/n : j = 1, ... ,n}. Then 
it is easy to check that the physical braid b of the motion of E has a mathematical braid 
(a-ia-2 . . . d-n-i)'^ and further that the diffeomorphism ^{z, t) — {\z\ exp(— 27rii), makes b 
weakly equivalent to htriv It is known that this mathematical braid and its powers are the 
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only cases where this happens ([Bml]). More precisely, if we let be with the added 
relation that (o"iO"2 . . . (Tn-l)^ is the identity element, then the braid word of b is equivalent 
to the identity clement of B^ exactly when b is weakly equivalent to hfj-j.,.. Note that the 
center of B^ (i.e. all the elements which commute with every other element) is exactly all the 
powers of ((7i(72 . . . an-i)^, so one could more succinctly define B^ as B^ mod its center. 

For n > 2, Bn has infinitely many elements, and so there are many diffeomorphisms 
that are not isotopic to the identity; the simplest example is when E is three points and 
F is a rel E diffeomorphism with mathematical braid af. When diffeomorphisms are not 
isotopic to the identity, isotopy can be described using the identity case as follows. Two rel 
E diffeomorphisms / and g are isotopic if there is third diffeomorphism h, with h isotopic to 
the identity, and g — hof. If / and g are time one maps of fluid motions, this says that they 
are isotopic if we can accomplish the same advection as g by first allowing the advection for 
/ and then following it by a fluid motion that keeps the points of E flxcd. 

The connection of braids to this more general notion of isotopy requires the general def- 
inition of weakly equivalent physical braids. Two physical braids b° and b^, both with end- 
points E are weakly equivalent if there is a level preserving diffeomorphism \E' : x [0, 1] — > 
X [0, 1] that takes b*^ to b^ where "^{z,t) = (■'tpf(z),t) for some family of diffeomorphisms 
ipt with ipo = ipi = id. Using the result in the identity case, we have that two rel E advection 
diffeomorphisms are isotopic if and only if their physical braids with endpoints E are weakly 
equivalent which happens if and only if the corresponding braid words are equivalent in B^. 
In different language, the collection of rel E isotopy classes forms a group under composition; 
this group is isomorphic to Bn, where n is the number of elements in E. 

It is worth noting that when two physical braids with endpoints E are equivalent in the 
sense defined in §4.1, they are weakly equivalent (hence the use of the modifier "weakly"). 
This is easily seen using the corresponding mathematical braids. The equivalent physical 
braids have the same braid word in Bn, and hence in Bn, and thus they are weakly equivalent. 
The converse need not be true, the simplest example being the trivial braid and a physical 
braid having braid word ((Jict2 . . . cxn-i)". These represent the same word in Bn, but not 
in Bn, and hence the physical braids are weakly equivalent but not equivalent. 

The last step is to extend these ideas to allow different distinguished sets E so that vortex 
motions with different initial configurations can be compared. A rel E diffeomorphism / and 
a rel E' diffeomorphism /' are isotopic up to conjugacy if there is a homeomorphism h with 
h{E') = E and hf'h~^ is isotopic to / rel E. Thus / and g are isotopic after a change of 
coordinates. Two physical braids b and b' with endpoints E and E' are weakly equivalent if 
there exists a diffeomorphism h : ^ with h{ai{0)) = a'-(0) for all i (and so h{E) = E') 
and the physical braid {{h{ai{t)),t)} is weakly equivalent to b'. We then have that two 
advection diffeomorphisms are isotopic up to conjugacy if and only if their physical braids 
are weakly equivalent, which happens if and only if their mathematical braid words are 
conjugate in Bn- 

§5.2 The Thurston-Nielsen representative. The Thurston-Nielsen theory provides tools 

to analyse homeomorphisms using their isotopy classes. The theory constructs in each isotopy 
class a special homeomorphism which is now called the Thurston-Nielsen (TN) representa- 
tive. This TN-homeomorphism is in a precise sense the simplest map in the isotopy class. 
Simplest here means in both the topological and dynamical sense. Thus once we under- 
stand the topology and dynamics of the TN-homeomorphism in an isotopy class, we know 
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dynamical and topological complexity that must be present in every homeomorphism in the 
class. This subsection and the next contain an introduction to the Thurston-Nielsen theory 
that is targeted for our applications. In particular, only homeomorphisms of the punctured 
plane are considered here, or equivalently, homeomorphisms of the plane rel a finite set E. 
The Thurston-Nielsen theory holds for the much wider class of surfaces of negative Euler 
characteristic. The reader is urged to consult [T], [FLP] or [CB] for a more balanced and 
complete treatment. For a survey of the dynamical applications of the Thurston-Nielsen 
theory, see [Bd]. 

The TN-homeomorphisms are of three basic types. The first type consists of dynamically 
very simple maps called finite order (fo). These are defined by the property that for some 
m > the map composed with itself m times equals the identity map. The second type of 
TN-homeomorphisms are dynamically very complicated and are called pseudoAnosov (pA) 
homeomorphisms. These will be discussed in greater detail in the next subsection because the 
presence of a pA homeomorphism in an isotopy class has strong implications for the dynamics 
of an advection homeomorphism. The final type of TN-homeomorphisms are reducible, which 
means that there is a collection of simple disjoint loops with the property that the loops are 
permuted by the reducible homeomorphism. Cutting along the loops yields a collection of 
smaller surfaces and on each of these surfaces the reducible homeomorphism is either finite 
order or pseudoAnosov. 

A given isotopy class can only contain a TN-map of one type, and thus an isotopy class 
is called pseudoAnosov, finite order or reducible depending on what kind of TN-map it 
contains. Note that the conjugate of a TN-map is also a TN-map of the same type. Thus, 
using the results of the previous subsection, one may also speak of the TN-type of a braid or 
a braid type. There are many theorems which help decide the TN-type of a braid or isotopy 
class. Rather than catalog these results here it will suffice to note that there is a computer 
implementation of an algorithm due to Bestvina and Handel ([BH]) (cf. [FM] and [Ls]) 
which, given the braid word, decides the TN-type. In the pA case the algorithm outputs the 
isotopy invariant dynamical data described in the next subsection.-^ 

For planar regions the finite order braids and rel E isotopy classes are well known. If we 
let Rn denote rigid rotation of the plane by — 27r/n, then certainly K!^ = id, and so is 
finite order. We may consider Rn as the time one flow map of the motion given in complex 
coordinates by ht{z) — —\z\ex.-p(27rit/n). If E is the set of n*'* roots of unity, its physical 
braid under this motion is (3n '■= criO"2 . . . Cn-i- Note that f]^ is the identity element of 
Bn, and making this property hold is, in fact, the defining feature of Bn- Thus [3^ is a finite 
order braid and so is for any < m < n — 1. A classic result of Brouwer (see section 
8.2 in [Bd]) says that all finite order homeomorphisms of the plane are conjugate to rigid 
rotations. This implies that the are the only finite order braids on n strands. Examples 
of other braids with various TN-types are given with vortex examples in §6.2. 

For planar regions the Thurston-Nielsen theory only has content for regions having three 
or more punctures or holes, or equivalently, for isotopy classes rel three or more points. 
When E contains just two points, all TN-homeomorphisms are finite order since B2 only 
contains the identity element e and the braid ai with af = e. Classic results of Alexander 
(see Cahpetr XV in [D]) say that when E has one or zero elements, there is only one isotopy 

^ There is a C-|— I- implementation of this algorithm by Toby Hall (with a Win95 interface) 
available for download. Contact the first author for site locations. 
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class; all orientation preserving homeomorphisms of the plane are isotopic, as are those of 
the once punctured plane. 

§5.3 PseudoAnosov maps and isotopy classes. Linear Anosov maps on tori provide 
one of the best known and understood examples of chaotic dynamics. The homeomorphism 
(j) of the two torus induced by the matrix 



acting on the plane is called Thom's toral automorphism or Arnolds' cat map, the latter 
nomenclature derived from an illustration in [AA]. The homeomorphism has a Markov 
partition with transition matrix M that gives a symbolic description of the dynamics. The 
largest eigenvalue A > 1 of M controls much of the dynamical asymptotics. The num- 
ber of fixed points of 0" grows like A" as does the length of topologically nontrivial loops 
under iteration. In addition, these properties are shared in the appropriate sense by any 
homeomorphism on the torus that is isotopic to ([Fk])- 

For the map there is uniform stretching and contraction by A at every point, and 
the corresponding eigendirections fit together to yield a pair of foliations of the torus by 
unstable and stable manifolds. Such nonsingular foliations cannot live on surfaces of negative 
Euler characteristic like the multi-punctured plane, and thus these surfaces cannot support 
Anosov maps. However, pseudoAnosov maps do exist on these surfaces and share many of 
the properties of Anosov maps. 

In pseudoAnosov maps there is still uniform stretching and contracting at each point by 
a factor A. The unstable and stable directions still fit together into invariant unstable and 
stable foliations but these foliations must have a finite number of points (called singularities 
or prongs) where there is either one or three or more stable directions and the same number 
of unstable directions. Despite the presence of these singularities, pseudoAnosov maps share 
most of the basic properties of Anosov maps. There is a Markov partition with transition 
matrix M that allows one to code the dynamics of the pA map 0. The largest eigenvalue 
of M is the stretching factor A > 1. The number of fixed points of 0"^ grows like A"^ as 
does the length of topologically nontrivial loops under iteration. The precise way in which 
these properties are shared by any isotopic map is described by the following theorem due 
to Handel ([H]). The notation g^y means the map g restricted to the set Y. 

Theorem: If cf) : — > is a pseudoAnosov homeomorphism and g is isotopic to 0, 

then there exists a closed, g invariant set Y , and a continuous, onto map a. : Y ^ M^, so 
that a o = o a. Further, for any periodic point x of 0, there is a periodic point of the 
same period y of g with Q:{y) = x. 

Thus any map g isotopic to a pA map has a compact invariant set that is semi conjugate 
to the pA map. Since the map a is onto, the dynamics of g are at least as complicated as 
those of the pA map 0. Thus, for example, the exponential growth rate of the number of 
fixed points of g'^ is at least A and the topological entropy of g is greater than or equal to 
log(A). Note that nothing prevents the dynamics of g from being more complicated than 
that of 0, the theory provides just a lower bound. 

A loop in a surface is called topologically nontrivial if it cannot be deformed into a point, 
a puncture, or a boundary component. Another basic result about pA maps is that the length 
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of topologically nontrivial loops grows like A" under iteration. Via the semiconjugacy, this 
implies a similar growth under g. It is important to note that this does not imply that the 
map g has Lyaponov exponents equal to log(A) everywhere. One does obtain from smooth 
ergodic theory that g has an ergodic invariant measure with an exponent at least log(A) (cf. 
[KH]), but this measure will be supported on the set Y, which could be small with respect 
to the usual measure on M. 

Since the number A can be computed by the Best vina- Handel algorithm, the braid de- 
scribing an isotopy class yields quantitative information about any homeomorphism in the 
isotopy class. An example of an isotopy class containing an advection homeomorphism in- 
duced by vortex motion is given §6.3. The Bestvina-Handel algorithm also returns a one 
dimensional graph called a train track in addition to a self-transformation of the graph. The 
edges of the graph form the Markov partition for the pA map, and the self-transformation of 
the graph provides the transition matrix as well as the structure of the invariant foliations. 
These foliations arc sometimes called the "invariant manifold template" of the pA map. The 
semiconjugacy then gives a lower bound or skeleton for the invariant manifold templates of 
9- 

§6 Advection and isotopy classes. 

In this section we use the Thurston-Nielsen theory to study the of advection of a passive 
particle in the velocity field produced by three point vortices. In §6.1 a Poincare map 
(or advection homeomorphism) is defined whose iterates describes the dynamics of passive 
scalars. This is a standard construction, but an adjustment is required because of the 
dynamic phase in the periodic motion of the vortices. In §6.2 we study isotopy classes of the 
advection homeomorphisms in the plane and array case using the braids computed in §4. 
Since all the motions in a single regime have the same braid, all the corresponding advection 
homeomorphisms have the same TN type. For the planar case, these are all reducible or 
finite order, but the array case has regimes of pA type. In §6.3 we study one such pA regime 
in detail. 

§6.1 The Poinceire map. As we have seen, a vortex motion Za{t) can be obtained by 
solving (2.1) directly, or as in §3, by using (3.4) and a period P solution Z{t) of the reduced 
equations (3.5). The vortex motion, in turn, generates a velocity field given by the right hand 
side of (2.3). The goal of this subsection is to define a Poincare map. The iterates of this map 
will describe the time evolution of passive particles advected in this velocity field. There is a 
standard construction of such a map for an equation of the form z — F{z, t) with F a velocity 
field P-periodic in t. The time-P fiow map / as defined in §5.1 by f{zQ) = z{P;zq,0) will 
satisfy /"'{zq) = z{nP; zq, 0), and so the iterates of / describe the time evolution of points 
in discrete steps of size P. 

For the advection caused by the point vortices this construction needs to be altered 
slightly because the vortex trajectories satisfy Za{t + P) — Za{t) + 6, with b the dynamic 
phase associated with Z{t). The easiest alteration is to pass to a uniformly translating frame 
with velocity b/P. If F{z,t) denotes the right hand side of (2.3) corresponding to the given 
vortex motion, define the advection iiomeomorpliism f as the time P-fiow map of the equation 

C^Fic + ^t,t)-^. 

This advection homeomorphism satisfies f^{zo) = z{nP;zo,0) — nb, and so / describes 
the position of a particle after advecting for one period P and then translating by —b. 
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Alternatively, this could have been taken as the definition of /, i.e. /{zq) = z{P; zq.O) — b 
with z the solution of (2.3). Points that are periodic under / correspond to advected particles 
that are translating along with the vortices, but after some integer multiple of P they have 
returned to their initial positions relative to the vortices. Although the velocity field is not 
defined at the positions of the vortices, the homeomorphism / can be extended (continuously, 
not differentially) to make the three initial positions of the vortices, Za{0), fixed points. 

In §4.3 the motion of a vortex array in the cylinder was transformed to the plane minus 
the origin by the conformal map T. The advection homeomorphism in this case is trans- 
formed to / = TfT~^. Note that / may be extended to make a fixed point. Thus all the 
Sa{0) from (4.4) can be considered as fixed points of /. 

§6.2 The TN type of regimes. As seen in §4.2 and §4.4, all the periodic orbits in the same 
regime in the reduced plane give rise to the same braid. By results in §4.1 this implies that 
all the corresponding advection homeomorphisms have the same braid type, which allows us 
to speak of the TN type of a regime. If this type is pA, then as in §5.3 all the advection 
homeomorphisms will share a certain set of dynamics. 

In the case of three vortices in the plane, only two braid types arose as descriptions of 
regimes. In one of them, a pair of vortices rotate around one other while the other vortex 
is uninvolved. This is a reducible class, with a finite order component containing the pair of 
interacting vortices. The other regime, D, had braid {(Ji(T2)^, which is the identity element 
in ^3, and the advection homeomorphism in this case is isotopic to the identity, with the 
isotopy consisting of unwrapping the plane via single rotation; this topologically undoes 
whatever advection was caused by the vortices. Since only pA classes or components yield 
dynamical information from the Thurston-Nielsen theory, the theory gives no information 
about advection dynamics in the planar case. 

The Thurston-Nielsen theory does provide dynamical information about advection in- 
duced by vortex arrays. The Best vina- Handel algorithm was used to compute the TN type 
of the various regimes from Figure 3.2. Using the braids given in the second column as input, 
the computed TN type is shown in the third column of Table 4.1. When a regime is pA (or 
has a pA component), the expansion constant is given in the right column. 

As we comment on the various regimes the reader is urged to draw a picture of the 
appropriate braids. The regimes which contain poles. III through VIII, correspond to a pair 
of vortices which rotate around one other while the other is uninvolved. As in the planar 
case, the isotopy classes of the advection homeomorphisms are reducible with a finite order 
component containing the pair of interacting vortices. 

In the braid for regime X shown in Figure 6.1a, the three right most strands do not 
link with the left most string. This shows that the corresponding isotopy class is reducible 
with a reducing circle containing the three right most strands. Inside this reducing circle the 
behaviour as a three braid is cr2cricr2cricr2cri = {aia2)^ using the relations in i?3. This braid 
is the same as the braid for region D in the planar case, which represents the identity class. 
Thus for regime X the isotopy classes are reducible with all finite order components. The 
braid for regime IX is reducible in a way similar to that of regime X, but inside the reducing 
circle the class is again reducible, but still with all finite order components. 

In the braid for regime II shown in Figure 4.4d. The schematic of the plane motion 
in Figure 4.4b makes it clear that the class is reducible; one can enclose the motions in 
successive, nested reducing curves. The classes for regime I are similar. 
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Figure 6.1: Mathematical braids of regimes from Figure 3.2: (a) 
Regime X (b) Regime XII (c) Regime XI. The braid of regime XI 
restricted to vortices 2, 3, and 4 is shown in (d), and a schematic of 
the plane motion of the vortices is shown in (e) 

The braid for regime XII is shown in Figure 6.1b. The braid makes it clear that vortices 
1 and 3 have essentially the same motion with respect to vortex 2 and the origin ("vortex" 
4). They do rotate around each other once, but they may be enclosed in a reducing circle. 
To analyse the class outside the reducing circle we think of squeezing the strands of vortices 
1 and 3 into a single strand, or equivalently, we delete one of the strands. The result is the 
three braid (tIo-iCtIo'^ shown in Figure 6. Id. This braid corresponds to a pA class rel 3 points 
with A ^ 13.9. 

Regime XI also yields a pA class with A fa 13.9. Although this class is irreducible, the 

strand corresponding to vortex 1 is in fact a fixed point of the pA map in the class defined 
by the other three strands. The latter is the same as that for regime XII (see Figure 6.1). 
The dynamics in regime XI are considered in detail in the next subsection. The last regime 
is XIII which is pA with A fa 9.9. 

§6.3 An example of the dynamical implications of the Thurston-Nielsen theory. 

As an example we describe some of the results afforded by the Thurston-Nielsen theory 
for advection homeomorphisms coming from vortex motions in regime XI. Recall that this 
regime's TN-type is pA with A ~ 13.9. By fixing a periodic orbit Z{t) in this regime and 
an initial position zi{0) we obtain a motion of all the vortices as in §6.1. The corresponding 
advection homeomorphism is denoted /. We may variously consider / as a homeomorphism 
of the singly-periodic plane, the cylinder, or as transformed from the cylinder to the plane 
by the conformal map T. 
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Using the theorem of §5.3 and the comments after it we know that the number of fixed 
points of grows with n at least hke 13.9"". These fixed points for f^^ correspond in the lab 
frame to passive particles that are advected by the vortex velocity field and after n periods 
have returned to the same position relative to the vortices. They have thus translated by n6, 
where h is the relevant dynamic phase. We also know that under advection, topologically 
nontrivial curves grow in length at least at a rate of 13.9"^. In addition, the topological 
entropy of / is at least log(13.9). This implies that / has a hyperbolic invariant measure 
(in the sense of [KH]) with Lyapunov exponents log(13.9), but no information is available 
regarding its support. 

As noted at the end of §4.2, although all the advection homeomorphisms arising from 
vortex motions in regime XI will have these properties, the time scale of the advection 
described by the homeomorphism will vary greatly amongst motions in the same regime. In 
addition, variations in the dynamic phase will make a significant difference in the motion as 
observed in the lab frame. 

The picture of the trajectories of the vortices transformed into the plane is also instruc- 
tive. Figure 6.1c shows the mathematical braid of this regime. Figure 6. Id shows just the 
motion of vortices 3 and 2 and the origin (= "vortex" 4), while Figure 6.1e shows a schematic 
of this motion in the frame of the vortex 2 after it has been transformed to the plane. The 
feature of note is that the trajectory of vortex 3 crosses itself in an essential way. We cannot 
remove this self-intersection by deforming the path in the complement of vortex 2 and the 
origin. This essential self intersection in fact implies that the isotopy class of the advection 
homeomorphism generated by vortex 3, vortex 2 and the origin is pA. The self- intersection 
forces regions of the fluid to be repeatedly pushed across pieces of themselves as they are 
advected. This forces the advection homeomorphism to constantly stretch and fold in an es- 
sential topological fashion. It is well known that stretching and folding is a basic mechanism 
for the generation of chaos and in this case it is topological in nature and is present in any 
map in the isotopy class. 

§7 Summary and conclusions. 

Three-vortex systems with zero net circulation have shown themselves to be worthy of 
both physical and mathematical interest. Using a combination of dynamical, geometric and 
topological methods a fairly complete analysis of these systems has been achieved. After 
fixing a value of the integrals, the reduction process maps the system onto a one degree 
of freedom Hamiltonian. This Hamiltonian system has a natural subdivision into regimes 
and we have seen that all the solutions corresponding to a single regime yield topologically 
identical vortex motions as described by their braids. Viewing the vortices as the "large scale 
motion" the use of braids allows one to make precise sense of the topology and structure of the 
large scales. In addition, the notion of reducibility of a braid or isotopy class allows precise 
description of the hierarchy within the large scale motions. For certain regimes the Thurston- 
Nielsen theory uses the topology of the vortex motion to obtain very strong conclusions about 
the advection caused by the large scales. Although point vortices describe a very idealised 
fiuid, it certainly of value to have model systems for which commonly used phrases such as 
"structure" and "large scale" have a precise meaning with computable implications. 

Although this paper has focused on a particular kind of vortex system it should be clear 
that the methods employed have much wider applicability. Any system involving the periodic 
motions of points in a surface has a natural description using braids. If there is a natural 



31 



homeomorphism on the complement of the points, then it can be studied using Thurston- 
Nielsen theory. The prototypical such system is an iterated homeomorphism of a surface 
in which case the collection of points is a collection of periodic orbits. This application is 
fairly well studied, see [Bd] for a survey. A closely related example is a periodically forced 
oscillator, see [McT] for a survey. Within fluid mechanics, a two dimensional fluid region 
with a periodic stirring protocol by rods is studied from this point of view in [BAS]. 

It is also clear that the methods of this paper apply to a wider class of vortex systems. 
For example, the methods presented here immediately apply to the periodic motions of n > 3 
vortices. However, in the n > 3 case, the reduced dynamics is described by a n — 2 > 1 degree 
of freedom Hamiltonian system. In this case one expects that the generic bounded orbit is 
not periodic, but rather only recurrent (via the Poincare Recurrence Theorem). At a time of 
near periodicity, such an orbit can be closed in a canonical way and its braid computed. The 
corresponding advection homeomorphism may then be studied using the Thurston- Nielsen 
theory. The theoretical task is then to understand how the conclusions derived from the 
closed motion can be applied to the advection caused by actual recurrent motions. Similar 
methods coTild be then applied to the motions of regions of very concentrated vorticity in 
two-dimensional turbulence and the resulting fluid mixing. 

A distinct advantage of the application of the Thurston- Nielsen theory is that it gives 
quantitative, a priori lower bounds for the chaotic dynamics based only on qualitative, 
topological data. However, the methods are only applicable after a punctured plane is 
obtained, and then the homeomorphism under question must be in a pseudoAnosov isotopy 
class. When the methods are apphcable an important question is the relative importance of 
the dynamics dictated by the Thurston- Nielsen theory compared to other sources of chaos. 
For example, in advection due to point vortices, the singular vorticity gives rise to very 
violent local stretching. As the vortices move this highly stretched fluid can fold over itself 
in the classic "chaos engine" scenario. While this effect can be deformed away by an isotopy 
relative to the vortices, preliminary numerical investigations indicate that in certain cases 
this effect is an order of magnitude greater than the lower bounds given by Thurston-Nielsen 
theory. For a smooth distribution of vorticity, this local stretching effect would diminish, 
and so the global effects given by Thurston-Nielsen theory would presumably gain more 
importance. 

An example where the Thurston-Nielsen effects are clearly dominant is the slow stirring 
of a viscous fluid studied in [BAS]. Three stirring rods in glycerine were moved in a protocol 
corresponding to pseudoAnosov braid. Dye tracers indicate an advection dynamics very 
close to the pA model. The structure of the invariant manifold template of the pA map is 
clearly visible. The fluid motion in this case seems to realize the least dynamical complexity 
allowed in its isotopy class which corresponds to very efficient stirring. 

The topological methods employed here are very successful in illuminating certain fea- 
tures of the vortex system. As with any method, there are situations of greater and lesser 
applicability, but when applied with care these topological techniques provide a powerful 
tool in the study of two-dimensional chaotic dynamics. 
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